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The system under consideration is a two-dimensional one-component plasma in fluid 
regime, at density n and at arbitrary coupling F = j3e^ (e = unit charge, j3 = inverse 
temperature). The Helmholtz free energy of the model, as the generating functional for 
the direct pair correlation c, is treated in terms of a convergent renormalized Mayer di- 
agrammatic expansion in density. Using speciflc topological transformations within the 
bond-renormalized Mayer expansion we prove that the nonzero contributions to the regu- 
lar part of the Fourier component of c up to the fc^-term originate exclusively from the ring 
diagrams (unable to undertake the bond-renormalization procedure) of the Helmholtz free 
energy. In particular, c(fc) = — F//c^ + F/(87rn) — /c^/[96(7r?T,)^] + 0{k'^). This result flxes 
via the Ornstein-Zernike relation, besides the well-known zeroth-, second- and fourth- 
moment sum rules, the new sixth-moment condition for the truncated pair correlation h, 
n{TrTn/2f J r^h{r)dr = 3{T - 6)(8 - 3F)/4. 
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1. INTRODUCTION 

Coulomb plasmas are the model systems for studying the effect of long-range interparticle 
interactions on statistics of classical lattice and continuous fluids. It was observed that, 
in arbitrary dimension, the long-range tail of the Coulomb potential gives rise to exact 
constraints, sum rules, for truncated particle correlations (for an exhausting review, see 
Ref. 1), namely the zeroth- and second-moment conditions^^'^'' . 

The concentration on two dimensions (2d) with logarithmic interparticle interactions 
and on the one-component plasma (OCP), i.e. the continuous system of charged parti- 
cles embedded in a spatially uniform background, brings some physical peculiarities and 
relevant mathematical simpliflcations providing additional exact information about the 
system, like: 

- A formal relationship to the fractional quantum Hall effect'-^^; 

- An experimental evidence for the Wigner crystallization at low temperatures'-^^ ; 

- The dependence of the statistics on the only parameter - coupling constant F ~ 
1 /temperature (the charge density scales appropriately the distance); 

- The availability of the equation of state'-^-* ; 

- The mapping to free fermions at special coupling F = 2^^'' (for various sample's ge- 
ometries, see review 8) characterized by a pure Gaussian decay of the truncated pair 
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correlation h; the evaluation of the leading term of the (F — 2)-expansion of h, indicating 
the change from monotonic to oscillatory behavior just at F = 2^'''); 

- The rigorous derivation of the weak-couphng F — > Debye-Hiickel hmit^^^ , and the sys- 
tematic F-expansion of h in terms of a renormalized Mayer diagrammatic expansion'^^'^'^^^; 

- The shift of the compressibility equation to the fourth-moment condition^^^"^^^ under 
the assumption of cluster conditions; 

- A symmetry of thermodynamic quantities with respect to a complex transformation of 
particle coordinates^^^^ for arbitrary coupling F, implying a functional relation among the 
pair correlations. The last is equivalent to an infinite sequence of sum rules relating tlio 
coefficient of the short-distance expansion of two-particle correlations (the lowest level 
of the sequence was derived by Jancovici*^^^''). The generalization of the symmetry to 
multi-particle densities, possessing specific invariant structure, was given in Ref. 17. 

- The suggestion that, at arbitrary F, the 2d OCP is in the critical state^^^'^^^ in terms of 
the induced electrical-field correlations (but not the particle correlations). The free energy 
is therefore supposed to exhibit finite-size correction predicted by the conformal-invariance 
theory, as was verified by the rigorous finite-size treatment of the F = 2 case and also 
numerically(^°^ , by using exact finite-size techniques^^'^^"^^^ , for coupling strengths F = 4 
and F = 6. 

The present paper is devoted to a rigorous derivation of the new sixth-moment sum 
rule for the truncated pair correlation h of the 2d OCP. The mathematical basis comes 
from the convergent bond-renormalized Mayer expansion in density^^^-*. Within a specific 
classification of the diagrams in the renormalized format for the Helmholtz free energy, 
the functional generator for the direct pair correlation c, we prove that the regular part of 
the Fourier component of c is determined up to the fc^-term solely by the (unrenormalizcd) 
ring diagrams of the generating free energy. This result implies via the Ornstein-Zernike 
(OZ) relation, besides the known zeroth-, second- and fourth-moment*^^^"^^-* sum rules, 
the explicit formula for the sixth moment of h. 

The paper is outhned as follows: 

In section 2, we recapitulate briefiy the ordinary Mayer diagrammatic representation 
in density for h and c pair correlations and the (excess) Helmholtz free energy as the 
generating functional for both of them. 

Section 3 deals with exactly solvable cases or limits of the 2d OCP, expressed in terms 
of /i-moments. These involve the momentum sum rules, the F = 2 coupling together 
with the leading (F — 2) correction term and the Debye-Hiickel F — >^ limit. Here, we 
have registered a very important fact. The three known (appropriately rescaled) zeroth, 
second and fourth moments turn out to be F-polynomials of (successively increasing) 
finite order. By using the renormalized Mayer expansion*^^^) we were able to compute 
for the as-yet-unknown sixth /i-moment the coefficients to a few lower orders of the F- 
expansion terms around the Debye-Hiickel limit and observed, within the range of these 
low orders, the finite F-truncation also for this case. To our surprise, this finite truncation 
represented an exact interpolation between the F ^ and F = 2 couplings and, moreover, 
reproduced correctly the leading (F — 2) correction term^^-*. Since the coupling F = 2 and 
its neighborhood do not play any special role in view of the F-expansion, the above fact 
was a strong indication and motivation for us. 

Section 4 describes the formalism of the renormalized Mayer expansion*^^^^ The renor- 
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malization of bond factors consists in a multiple-bond expansion of Mayer functions and a 
consequent series elimination of field circles, resulting in the modified Bessel functions of 
second kind. The novelty lies in the classification of diagrams representing the Helmholtz 
free energy - the functional generator of c - according to the possibility of performing the 
series-elimination transformation: (1) simple unrenormalized bond generates the charac- 
teristic singular term of c; (2) all unrenormalized ring diagrams (which cannot undertake 
the series-elimination procedure) generate the renormalized "watermelon" Meeron graph 
contribution to the regular part of c; (3) every other diagram is expressible with all bonds 
renormalized and as generator gives rise to a family of c-diagrams: the families do not 
overlap with one another and, as units, they exhibit remarkable "cancellation properties" . 

The "cancellation" phenomena is the subject of the crucial section 5 where we prove 
that, regardless of the topology of a separate graph belonging to class (3), the zeroth 
and second real-space moments of the c-diagrams family, generated from the underlying 
graph, vanish. The proof of the nullity of the zeroth moment follows from a trivial scaling 
property of the Bessel functions with respect to density. The proof of the second-moment 
condition is much more complicated. Besides the above scaling property it requires to in- 
troduce an elimination procedure for two-coordinated root points generated on renormal- 
ized bonds and to reveal "hidden zeros" due to the translational and rotational invariance 
of the infinite system, realized through per-partes integration of field-point coordinates. 

In Conclusion, after evaluating the only contribution to the Fourier component of c 
up to the /c^-term, namely that of the renormalized Meeron graph, we write down by 
using the OZ relation the explicit formula for the sixth moment of h. The structure of 
higher-order moments of h is also discussed. 



2. A SKETCH OF THE MAYER EXPANSION IN DENSITY^^^) 

We consider a system of identical pointlike particles in volume V of a (i-dimensional 
space, interacting through pair potential v; v will occur in combination called the Mayer 
function 

/(i,j)=exp[-/3T;(i,i)]-l (1) 

where (3 is the inverse temperature and, for notational convenience, a position vector 
Tj is represented simply by i. In the inverse (density) format, i.e. with the density 
n(r) = (I]j5(r— rj)) as controlling variable, the excess free energy F'^^ is the relevant 
thermodynamic potential. Its Mayer diagrammatic representation reads 

ppex _ connected diagrams which consist of N >2 field (black) 

n— circles and /—bonds, and are free of connecting circlesj (2) 

(the removal of a connecting circle disconnects the diagram). The excess free energy is 
the generating functional for the truncated pair correlation 

n.(1.2)-n(lK2) 
^ ' ' n(l)n(2) ^ ^ 
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with the two-body density n2(r, /) = (I^j^j S{r— ri)6{r' — r,)) and for the direct corre- 
lation function c, in the sense that 

With regard to (H), this imphes 

h{l,2)= {all connected 1,2 — rooted diagrams which consist of field 

circles and /—bonds, and are free of articulation circlesj (6) 

[the removal of an articulation circle disconnects the diagram into two or more compo- 
nents, of which at least one contains no root (white) circle]; 

c(l, 2) = |all connected 1, 2 — rooted diagrams which consist of field 

n— circles and /—bonds, and are free of connecting circlesj (7) 

The link between h and c is established in terms of the OZ relation 

h{l, 2) = c(l, 2) + J c(l, 3)n(3)/i(3, 2)d3 (8) 

If the system is infinite (V oo), homogeneous, n(l) = n, and both isotropic and 
translationally invariant, h{l, 2) = h{\l — 2\) , c(l, 2) = c(|l — 2|), it is useful to introduce 
the Fourier components 

/(r) = \ f exp(ifc- r)/(fc)dfc (9a) 



(27r)'^/2 

Especially, in d = 2 dimensions, 

f{k) = / rf{r)Jo{kr)dr 



{2ttY/2 

/>) = 7^J^M-^k-r)f{r)dr (9b) 



f:7-^f?yi/-'V(r)dr (10) 



j=0 



ijiy \ AJ 2n 



with Jq being the ordinary Bessel function. In the Fourier space, the OZ relation (|^) takes 
the form 

h{k) = c{k) + {2T:Y'^nc{k)h{k) (8') 

where k = \k\. 
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3. MOTIVATION 

The classical OCP is a system of particles of charge e embedded in a spatially uniform 
neutralizing background. In = 2 dimensions, the Coulomb interaction energy is given 
by 

-pv{l,2) = rin|l-2| (11a) 

-(3v{k) = -r/e (lib) 

with r = /5e^ being the coupling constant. We will concentrate on the thermodynamic 
limit of the fluid regime with constant density n{i) = n and use the notation 

I^. = f r^ih{r)dr (12) 



for the moments of the truncated pair correlation. 

Let us first summarize exactly solvable cases of the model. In the weak coupling F — 
limit, h displays the Debye-Hiickel screening'^^^ 



h{r; r ^ 0) ~ -rKo{rV27crn) (13) 
with Kq the modified Bessel function of second kind. Consequently, 

lhnn(^)'j2,(r) = -(j!)^ (14) 

At r = 2, the mapping onto free fermions provides a pure Gaussian form of h^'^\ 

/i(r; r = 2) = — exp(— Tmr^) (15) 

implying 

n{nnyh,iT = 2) = -j\ (16) 
The leading-order of the series expansion around F = 2 results in^^^ 

h{r;T) = h{r;T = 2) + {T - 2)6h{r) + . . . (17a) 
6h{r) = Ei{~7inr^) - ^Ei{-7mr^/2) 

+ exp(-7mr2) |^Ei(7™rV2) - [\n{7mr^) + C]| (17b) 

where C is Euler's constant and Ei the exponential-integral function. From (|l3) one gets 
after some algebra 



n 



/2.(r) = -j\ + (F - 2)ji f E ^ - 1 1 + 0[(F - 2) 



The long-range tail of the Coulomb potential gives rise to exact constraints (sum rules) 
for the moments of the truncated two-body correlation, like the zeroth-moment (perfect 
screening) condition 

n/o = -1 (19a) 
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the second-moment (Stillinger-Lovett) condition*^^'^) 



/TcTn\ , , , 

n — - /2 = -1 (19b) 



the fourth-moment (compressibihty) condition'-^^ 



n[^-^ U = -^^,Y (19c) 



Note that the sum rules are consistent with exact formulae (14), (|TI 

The expressions for the rescaled moments (|19D correspond to finite truncations of their 
F-expansion around F = 0. The F-expansion technique of explained and extended in 
the next section, enables one to evaluate systematically the coefficients of the F-expansion 
also for higher moments. For the sixth moment of /i, we were able to attain with a little 
effort the third order of F, with the result 



n (^) ' h = -36 + |f - ^F^ + X F=^ + O(F^) (20) 

The appearance of the zero coefficient to the F^ power indicates the possibility of a finite 
F-truncation for Jg, too. This hypothesis is strongly supported by the fact that the 
truncation of (pO|) at the F^-term interpolates correctly from the F — > limit [relation 



(|T4D] to the F = 2 coupling in the sense that the rescaled moment Iq satisfies (0), i.e., 
the rhs of (^) acquires the exact value (—6) at F = 2 and exhibits the exact prefactor 
21/2 to the leading (F — 2) correction term. These facts were our primary motivation for 
proving rigorously the conjecture of the finite F-series-truncation of Iq. 

It turns out to be more convenient to formulate the above sum rules in terms of the 
small- expansion of the Fourier component of the direct correlation, c{k). Inserting (^^) 
together with the suggested truncation of the sixth moment at the F^-term ( pOD into the 
representation (jT^) for h{k), the OZ relation (8') implies the expected form 

= + ^ - ^J-T2 + (21) 

Note the characteristic singular leading term of c(fc), —jSv^k) = —T/k"^, succeeded by the 
regular /c^-series expansion part: its knowledge up to the k"^^ term determines h{k) up to 
the /c^*^-'"'"^) term, or equivalently, the real-space /i- moments ( JT^ up to l2{j+2)- 



4. RENORMALIZED MAYER EXPANSION 

The Mayer function / can be expanded in the inverse temperature f3 as follows 

/(I, 2) = -Pv{l, 2) + ^[-Pv{l, 2)]2 + l^[-Pv{l, 2)f + ... (22) 
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Graphically, 



o o = o o + a x> + a 50+.. . (22') 

1 2 1 2 l"~-"2 1"~'"2 

where the factor 1/ (number of interaction lines)! is automatically involved in each dia- 
gram. Let us perform the above /-decomposition within the diagrammatic representation 
(D of 




^ ^ *■< ^ ^ 4 ^ >^ ^ ^ 



f- -t f- -t ^ _ ^ ^ ^ ^ ^ ^ 

I l_l l_l l_ 

etc. If there are only one- or two-coordinated field circles in a graph (under coordination 
of a vertex we mean the number of bonds met at this vertex), we do nothing. If there 
are some three- or more-coordinated field circles in a graph, we can eliminate all two- 
coordinated field circles by a series transformation and arrive at a connected graph of field 
circles of coordination > 3, called skeleton. Grouping the diagrams which are reduced 
to the same skeleton after series elimination, the bonds connecting skeleton field circles 
become dressed according to 

i^'(l,2) = a---o + o----»----o+ a---«----»----o-|-... 

12 1 2 1 2 

= cy-^^^-^-y^ (23) 

1 2 

Equivalent ly, 

K{1, 2) = [-pu{l, 2)] + I [-pu{l, 3)]n{3)K{3, 2)d3 (23') 
or, in the case of an infinite homogeneous fluid, 

K{k) = [-(3u{k)] + {2nf/\[-Pu{k)]k{k) (23") 
The procedure of bond renormalization thus implies 

oo 

/3F-[n] =•.--- H, + Do[n]+ J2 Ds[n] (24a) 

s=l 
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where Dq is the sum of all unrenormalized ring diagrams (which cannot undertake the 
renormalization procedure) 



Dn 



> + 



A 

/ \ 



and 



t- 

X + ; 

^ 4- 



-t 

; + 

-4 



(24b) 



s=l 



^ = I all connected diagrams which consist of > 2 field 

n— circles of bond — coordination > 3 and multiple 
K — bonds, and are free of connecting circles | 



(24c) 



represents the set all remaining completely renormalized graphs. Under multiple K- 
bonds we mean the possibility of an arbitrary number of i^'-bonds between a couple of 
field circles. The order of numeration of D-diagrams in (24c) is irrelevant, let us say 



(25) 





and so on. The symbol Dg will reflect the notation of a given diagram and simultaneously 
its integral representation. 

Having classified the renormalized graphs of (3F^^, we proceed by considering the 
direct correlation c, defined by equation (5). From (|2^) one derives 

oo 

c(l, 2) =c> - - o + co(l, 2) + Y. c,(l, 2) (26a) 
1 2 s=i 

where co(l,2) = 5n{l)5n{2) can be easily shown to correspond to the renormalized 

"watermelon" Meeron graph 



-o(l'2)= <Z>^ =J,K\l,2) (26b) 

and Cs(l, 2) with s = 1, 2, . . . denotes the whole family of 1,2-rooted diagrams generated 
from Dg, 

To get explicitly a given family c^, one has to take into account the functional dependence 
of the dressed -ft^-bonds (|23|) on the density as well. Since with regard to (23') it holds 

^f|^ = /<-(1.3)A-(3.2) (27) 



the functional derivative of Ds with respect to the density field generates the root circle not 
only at field-circle positions, but also on i^-bonds, causing their "correct K-K division". 
For example, in the case of generator Di drawn in (|25|) , we obtain 

2 1 

r^^^^'^^^ pj-^O^ r^^^^^^ 
iVw-^ l\ / 2Vw-^ 

1 1 2 

, ^ + ^ (28) 

2 

It stands to reason that the coordination of field circles remains to be > 3 after the 
functional derivation, while the root 1,2-circles can be two-coordinated (just when being 
generated on a i^-bond). The {cs}^i diagram families evidently do not overlap with each 
other. 

The specialization to the infinite 2d OCP, with dimensionless interaction energy (jT!]), 
leads to the dressed i^'-bond (E3f) of the form 



K(k) = -— — (29a) 

r r 1 



Kir) = / — — exr)(ik ■ r)dk 



= -rKo{rV27irn) (29b) 

The series renormalization of the logarithmic interaction thus leads to the modified Bessel 
function of second kind: its decay to zero at asymptotically large distance makes the 
renormalized Mayer diagrams properly convergent'-^^-'. Note the specific r\/n dependence 
of K which has a fundamental impact on the n-classification of renormalized diagrams. 
As concerns the F-order of a given diagram D with field circles and L bonds, every 
dressed bond (p9D brings the factor F and enforces the substitution r' = r^/T which 
manifests itself as the F~^ factor for each field-circle integration ~ / rdr, so that the F- 
order = L — N. For example, in (p5|), Di ~ F and D2, D^, D4, constitute the complete 
set of /3F^^-diagrams ~ F^. 



5. SUM RULES 

Let the given completely renormalized diagram (s = 1, . . .) of the excess Helmholtz 
free energy be composed of skeleton vertices i = 1, . . . ,N and L bonds a = 1, . . . , L. 
Ds can be formally expressed as 

. N L 

Ds[n] = / n ^(0] n K{aua2) (30) 

•' i=l a=l 
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where ai, 02 G {1, . . . , A^}, ai < 02 denotes the ordered pair of vertices joint by the a- 
bond. Whenever not confusing, we will use the symbol Ka = K{ai,a2) and omit the 
ranges {1,A^} and {1,L} for a sum or product over skeleton vertices i and bonds a, 
respectively. The family of direct correlations Cs{r, r'), generated according to Cs{r, r') = 
S'^Ds[n]/Sn{r)6n{r^), can be straightforwardly expressed in the uniform- density regime 
n{i) = n as follows: 

c,(ry) = n''-'Jl[dzJ2 6{r-z)6{r^-j)Y[K^ (a) 

i i,j ct 

+n^-i / n E ^(^ - E a^) n (6) 

+n^-i / n E ^(^' - ^) E ^)K{r, a,) [] (c) 

/ X{diY.K{a,,7^)K{rJ,r)K{r,a2) W (d) 



/ n E ^)^(^' '^)^(^' «2) n (e) 

+n^|nd^ E i^(ai,r-)i^(r,a2)my)^(r',/52) 11 ^7 (/) 



where we have applied the functional relation (0). The (a) term on the rhs of (31) 
corresponds to the creation of root points at two skeleton vertices, the next two (6, c) 
terms to one root circle generated at the skeleton and the other one at a bond, the (c?, e) 
terms to two root points at the same bond and the last (/) term represents two root 
points generated at different renormalized bonds a ^ (3. It stands to reason that now 
K{r, r') = K{\r— r'\) satisfying (Al),(A2a,b) and, consequently, Cs(r, r') = Cs{\r — r^\). 
In this section, we aim at proving the validity of the moment equalities 



Cs{r)dr = (32a) 
J r^Cs{r)dr = (32b) 

for every family s = 1,2, . . ., regardless of the topology of the generating diagram Dg. 
To keep the interchange-particle symmetry and the translational-invariance property of 
the problem, we will use instead of (32a,b) the following equivalent definitions of the 
moments: 

# = ^Jcs{r,rJ)drdr^ (33a) 
4'^ = ^ J\r-r^\^Cs{ry)drdr^ (33b) 

The introduction of volume V into the formulation of the infinite- volume limit does not 
mean any loss of rigour. Definitions (33a,b) have to be understood in the sense that an 
arbitrary one of integration coordinates r, ■r', {i} can be taken, due to the invariance of 
the integrated function with respect to a uniform shift in all variables, as a reference and 
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put at the origin 0, with the simultaneous cancellation of V. The right choice of the 
reference can simplify otherwise tedious algebra. 



5.1. Proof of the zeroth- moment condition (32a) 

By using the definition (33a), the zeroth moment of (31) can be expressed as 



J, 



(a) 
(b + c) 
id) 
(e) 
(/) 



1+01,1} 



where we have taken into account relation (A2a). Since 
/dr 



K{r,a2) + K{ai,r) 



dn 



dn 



we find 



jis) _ ^ 



n 



N , N 



i=l a=l 



(34) 



(35) 



Let us put say i = 1 at the origin 0, and "cancel" the integration over 1 with volume V. As 
K{ai,a2) = -rXo(V27rrn| Oii — a2\) for the 2d OCP, the evoked substitution i' — WinVn 
for [N — 1) remaining coordinates i — 2, . . . ,N results in the factor l/n^^~^\ Therefore, 



jis) 
'0 



N 



- 



(36) 



in agreement with (32a). 
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5.2. Proof of the second-moment condition (32b) 

By using the definition (33b), the second moment of Cg (31) is written as 

N-2 I. 

/^d^EK-J^^^a 



X 



n 



V 

V 
2n^ 



J n^^E / \r-t\^K{ai,r)K{r,a2)drll Kf, 
f n^^E / \r- r'\'^K{ai,r)K{r,r')K{r',a2)drdr' [] K/s 

xK{P,, T')K{r', /32)dniT^ R^^a,^ K, 



V 



(a) 

{b + c) 
(d + e) 



a, 9 
("5^/3) 



(37) 

The integrations over rand r' in (37) correspond to root points generated on renormahzed 
X-bonds, and interacting with another root point at the skeleton or with one another. 

In Appendix A we show how to transform these integrals to the form with exclusively 
skeleton \i — jp,$(|i — j\) and ^(|i — j\) interactions [for definitions of $ and ^ see 
(A7) and (AlO), respectively] and appropriately "decorated" bonds [under decoration, we 
mean the derivation with respect to density, equation (A2a)]. The successive application 
of formulae (A8), (AlO) and (A9) to the respective terms (b+c), (d+e) and (f) of (37) 
yields 



n 



N-2 



V 

2n^-i 



+ 



V 
V 



Rd^E 



a L i 



E^(l«l-^r + K-^lO^ + ^^*« 



J /3^a 



Ild^E^all^/^ 



(a) 

(6 + c) 
{d + e) 



n f 1 

+y- yild^ E [^{\c^i-Pi\' + \c^i-/32\'+\a2-/3if + \a2-(32\') 



("5^/3) 



X 



dn dn dn ^ dn °'\ ^i- ^ 



if) 

(38) 



Let us now define the auxiliary function 

2n^+i 



n 



V 

N+l 



V 



(39) 



Up to the prefactor IjV ^ it originates from Dgin) (30) by picking out successively bonds 
one after the other and interchanging the bond factor K — > is expressible from 
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(A7) and (29a) in the form 

4r2 



(27rrn)2y (p2 + i)4 



exp 



ip ■ (cti — a2)\/'2,nTn dp 



(40) 



With the aid of the scahng argument leading to relation (36), G^*^ can be shown to scale 
with n like 



n 



N+l 



— n 



As a consequence, dG^^\n)/dn — 0. Exphcitly, 

2{N +l)n^ 



= 



+ 



V 



n 



V 

N+l 



V 



(«7^/3) 



Subtracting {Eq.(42)}/n from formula (38), the latter takes a simpler form 

,Ar-2 



j: 



is) 



n 



V 



n 



+ 



V 
2n^ 



V 



-/nd^E(i«i- 

i i,c 



^1 + |q;2 — "i 
ia(n$o 



2^ "^^^ n ^1 



dn 



n dn 



+^ / lid* E (l«i-A|' + |tti-/52r + |a2-/3ir+|a2-/92 



a, 13 



X- 



/3 



n 



dn dn 



It is trivial to show that 
1 a(n$«) 



_ — g*P-{"i-a2)^ 

n dn J 



TT 



= |Q;i — q;2| 



ip-{a\—a2) 



An{p^ + 27rrn)2 (p2 + 27rrn)3 

r2 7rr3 



dp 



4n(p2 + 27rrn)2 (p2 + 27rrn)3 



(41) 



(42) 



(43) 



dp 



4^2 



(43') 

Equation (43) together with the complementary one (43') represent the formulation of 
the second moment of Cg in the renormalized format in terms of quadratic interactions 
exclusively between pairs of skeleton vertices. 
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For tactical reasons, we now specify the connectivity of diagram D^, instead of enu- 
merating present renormalized bonds a = 1, . . . , L, by the set {z/jjl-^^^i where z/jj = Vji is 

the number of i^-bonds between skeleton vertices i,j {uij = if there is no bond between 
Let us choose a couple of skeleton vertices, say 1 and 2, and group all factors in 
(43), (43') associated with the |1 — 2p interaction: 



n 



N-2 



V 



/nd.ii-2r n K 



Ti ,■(; — 3 
{u<v) 



N N 

2K'^«(i,2)n^'^"(i,on^''^(2,j) 

i=3 j=3 



n 



dK{l,2) 
dn 



N 



N 



i=3 j=3 



n—\{K'''^{l,i 



dn 



i=3 



N 



N 



+ K'^\l,2)X{K'''^{l,i) 

1=3 



i=3 

Q N 



dn 



i=3 



1 



n 



dn 



+ ^Z/i2(z^l2-l)/^"^^-'(l,2) 



K{1,2 
dK{l,2) 



n 



dn 



N N 

ni^^"(i,^)n^'^'^(2,j) 

«=3 j=3 
2 Af N 

n^'^^'(i,on^''^(2,j) 

i=3 



j=3 



+ ^z/i2if^--i(l,2) 
+ ^z/i2i^^--i(l,2) 



n 



97^(1, 2^^ ^ 



Q iV 

i=3 



AT 



n- 



j=3 



n^''^(2,j 
^^n^"^(2,j) 

j=3 



+ 2^'^"(1,2) 



. ^"^ i=3 



9 ^ 

i=3 



(44) 



In the case of the considered Bessel functions (p9D , the operator n{d/dn) acting on K{i,j) 
can be substituted by a coordinate-operator as follows 



d 



dr. 



K{\n-r,\) 



^ -R,,K{\n-r,\) 



(45) 



As one can derive directly from the definition (^5|), there exist more equivalent represen- 
tations of operator Hij, 

R = r — 
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[rj - n) ■ V,- 



-(r,-r,)-(V,-V,: 



(46) 



Finally, denoting Fij = K°'^^{\i — j|) and summing over all pairs of skeleton vertices, Eqs. 
(43), (44) and (45) imply 



+ 
+ 



X 2 + 



i u<v i<:j 

— (Ri,-F„-) + 



-TT (R-ifc-^ifc) + -Tr~ (R-ifc-^jfc) 

-'^ ik ^ jk 



i^ikFik) + (RjkFjk) 

^ ik jk 



E 



k,l^i,j ^ikFjl 



[R-jfc-^jfc) (R-ji-^j 



(47) 



In what follows we aim at proving the nullity of the rhs of (^Tj), irrespective of the 
particular form of the bond-dependent functions Fijijij) (provided that the integrals exist 
what certainly applies to our case of F-functions) . As shown in Appendix B, due to a 
scaling transformation of coordinates in integrals of translationally-invariant functions 
over infinite 2d space, Eq.(|T7|) is reducible to a simpler relation 



X 



1 



E 



-— (n - Vj) ■ WiF.k + TTET^r-j - n) • WjFjk 

^J^ik ^J^jk 



+ E 



1 



k,l^i,j ^ikFjl 



[(r, - n) ■ V.F,fe] [(n - n) ■ V,F, 



(48) 



Let us consider {Fij} to be the functions of rfj rather than rjj, without changing the 
symbol F. Thus 



Using the notation 



Eq. 



can be rewritten as follows 



2\ 2 n - rj 



1 dF,,(4) 



F,j d{rt 



(49) 



(50) 



iFij 2(tj Vj^ Fij Fij 



(49') 
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Inserting (49') into (^Sp and grouping the F and FF terms we obtain 
8j(^)/n^-2 = /nd^ n Fu.[ E E ■ {rln, - r|,r,fc) 

i u<v (^i^j) ky^ij 

+ 2 E E FijFki [rU^j ■ rkj)iru ■ vm) + rUnj ■ rij){rik ■ nk) 



(i<j) (k<i) 



+r%{rji ■ rki)irji ■ Vki) + r^irji ■ ru){rjk ■ nt)] } (51) 



with the obvious notation = — r,-. 

Our further goal is to prove the nulhty of the rhs of (PI). We show that the per-partes 
integration of the bihnear FF term in (^) gives a hnear F term which exactly eliminates 
the linear F term in (|5l|). To do so, we express the bilinear FF term in integral ( ^l]) in 
an "ansatz" form 



n n Fuv \{ri2Fi2 + risFis + . . . + VinFin) ■ Qi 



+ (r2iFi2 + r23F23 + . . . + raAfFaTv) ■ Q2 + • • • 



(52) 



The vector functions must be linear combinations of Fki with k,l ^ i since the terms 

etc. In general, 



FijFij vanish in (|TD; = Qf'^^F^s + + • • • + 



l(2,4) 



JAf-l,Af) 



= E '^^^ 



kl 



(53) 



k<l 



where the three-point vector coefficients 

gf ''^ ^0 iff z ^ A; ^ / (54a) 
depend only on coordinates {rj, r^, r^}. It is natural to extend their definition as follows 

gf = g,f''^ (54b) 

Comparing Eq.(52) with (51), the coefficients to the Fi^Fu {i ^ k ^ I) and FijF^i {i 7^ 
j ^ k ^ I) terms imply the respective restrictions on g-vectors: 



fki ■ Qk Qi 



(55a) 



■ Qf''^ + r,, ■ gf + Vk, ■ Qt'" + nk ■ Q] 

rlki^ij ■ rkj)irii ■ r^) + rlinj ■ rij){rik ■ 

■ ni)irji ■ Vki) + rlirji ■ rii){rjk ■ nk] 



(55b) 



with no order inequalities put on {i, j, fc, /}. The important feature of the ansatz proposal 
(|52D consists in the equality 



n Kv {ri,iFa + ... + n^NFiN) = [ [] ^™ ) 

u<v \u<v / 



(56) 
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[see relation (49')]. Using formula (^) for every term in (|S2|), the consequent per-partes 
integrations lead to 

(52) = -^/nd^n^™EV.-Q. (57) 

The point is that the functions {F^i} in (|53|) do not depend on coordinate rj, and 
therefore the nabla operator acts only on coefficients { Qi'^'''^}- This is why expression (^) 
becomes linear in F-functions. It compensates the linear term in (^Tj) just when 

V,-Qf'''^ = 2ni-{rlni-rlru) (58) 

To summarize, J2*\ given by (0), equals to zero provided there exists a three-point 
vector function Q\'''^^ with properties (54a), (54b) and satisfying conditions (55a), (55b) 
and (|58D . One can readily verify on computer that such vector function exists: it is the 
homogeneous polynomial of the fifth order and its x-component reads 

2 

= 3 (-2M^fa; - 2m^w^ + 2ulvl + lulvl - 2u\vl - 2u\vl - u^v^ - u\v^ 

-uluyVy - Uyvlvy - 'iUn^vlvl " Sululv ^ " MajU^t^J, " UyV^Vy 

+6u^ulvl + Gulv^Vy - Qu^ulvl - Gulv^Vy + Su^UyVy + Sulv^Vy^j (59) 

where u = — Vi and v = ri — r^. The y-component [Qi''''^]y results from (^) under 
interchange transformation <-> Uy,Vx ^ Vy. We conclude that Jg*"* = 0, confirming 
relation (32b). 




6. CONCLUSION 

The Fourier transform of (26a) results in 

oo 

c{k) = -V/e + co(fc) + ^ c,(fc) (60) 

s=l 

With regard to the equalities (32a,b) proved in the previous section, we have 

cs{k)=0{k^) (s = l,2,...) (61) 

Using the formula*^^^-* 

x'^'^Kl{x)dx = ^^''-'^J^^_ is > 0) (62) 

the contribution of the renormalized Meeron-graph (26b) to the Fourier component of the 
direct correlation reads 

c,{k) = ^--^^+0{k') (63) 
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Consequently, the expansion of c{k) up to the fc^-term coincides with the suggested formula 
(pTD and, via the OZ relation, the (rescaled) sixth moment of h is indeed a finite F- 
truncation 

J r^h{r)dr=^(r -Q){8-3r) (64) 
as indicated in (pO|). 

In conclusion, we would like to stress that the derivation of the sixth moment of h was 
possible due to the property (pT]), valid separately for each family of direct-correlation dia- 
grams generated from the corresponding, completely renormalized, graph of the Helmholtz 
free energy. The higher-order coefficients of the fc-expansion of Cs{k) (s = 1, 2, . . .) we were 
able to attain do not longer vanish, e.g., 

gi(fc) = - (2^^)3 (4!)25! X K',{x)[176x+W8x' + 9x']dx + 0{k') (65) 

where the numerical value of the integral ~ 116.68 . . .. This fact indicates that the eight- 
and higher-order (appropriately rescaled) moments of h are probably infinite series in F. 
However, there might exist another mechanism for the exact solvability of higher-order 
moments as well. We believe that the present method will answer this interesting question. 
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APPENDIX A 

Here, we derive an algebraic procedure removing a two-coordinated root point, gener- 
ated by the functional derivation with respect to density at a renormalized -ft^-bond and 
interacting quadratically with another root point (generated either at a skeleton more- 
than- two-coordinated vertex or at a bond), in the representation of the second moment 
of a Cg-family (subsection 5.2). We consider the uniform regime with constant density n 
and translationally+rotationally invariant interactions 

f - dp 

K{ri, ra) = / exp[ip ■ (n - r2)]K{p)— (Al) 

J ZTT 

K as the function of n is given by the uniform analogue of functional relation (^Tf ) , 
or, in the 2d Fourier picture, 

-- 271 K\p) {A2b) 



dn 
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Let us first consider the case represented schematically as follows 



f{ui,U2)= ^ — ^ — = r'^K{\ui - r\)K{\r- uUiDdr {A3) 

M 

where Ui and 1*2 skeleton vectors define the bond decorated via { ^7\ ) by the two-coordinated 
root point r, integrated out and interacting via r^-interaction with another root point 
(whose coordination is irrelevant at this stage), put for simplicity at the origin 0. Using 
the Fourier representation of K, (A3) can be rewritten in the form 

= -i/dr/g/|/.(p)A-(,) j^^j 



Integrating twice per partes in (A4) over p (q) [what can be certainly done for our K (29a) 
with vanishing boundary contributions] and then integrating over r, implying S{p — q), 

we get 



(A5) 



(A6) 



fiu^,U2) = -y'dp|e^^("^-"^)K(p)Ai^(p) + ^(e-*''«^Ap[e**'-«^] 
+e^P-"iAp[e-*f "2]) + k{p)Vk{p) ■ Vpe'P("i-"2)| 

Integrating per partes once more the last term in (A5), we finally arrive at 

+ye*^'-("^-"^)|Vi^(p)|Mp 

With regard to the "decoration" relation (A2b), (A6) can be written in a more consistent 
form 

fiUl,U2) = -{Ul + UD ^^ + mm-W2\) ^^^^ 

$(w) = ye*P-"|Vi^(p)|Mp 
This equation admits a trivial generalization 

j \r— u\'^K{\ui — r\)K(\r— tt2|)dr 

= - {\u, - u\^ + \U2- u\^) ^ + $(|ui - val) 

The double application of formula (A8) solves immediately the problem of two root 
points generated on two different bonds: 



Vi r' V2 
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The final result reads 



J\r- r'\'^K{\ui - r\)K{\r- v^\)K{\vi - r'\)K{\r' - V2\)drdr' 



Ui - Vi\ + \Ui 

^K{\u,-U2\)^K {\v^-V2\) 
dn 



\U2- V2\ 



dn 



dn ^' ^ ^" ' dn 

When the two root points are generated on the same bond 



iA9) 



Ui r 



there holds 



K{\ui - r\)\r- r'\^K{\r- r'\)K{\r^ - 1*2! 



APPENDIX B 



In this part, we establish the transition from relation (^) to Eq.(|8). The origin of 
bond-dependent F-f actors is irrelevant with the only proviso that the integrals exist. 
Let us first consider the integral over infinite 2d space 



^ J In - ralVdn - r-2|)dridr2 = J r^f{r)d'i 



(51) 



written in the center-of-mass inertia; the supposed translational invariance of function 
fi2 = f{r 12) is important. The scaling transformation of coordinate r — >■ (1 + A)r does 
not alter the infinite boundary, and hence 



j r2/(r)dr= + j r"^ f[{l + X)r]dr 



Expanding (1 + A)^ = 1 + 4A + GA^ + 0{X^) and 



/(r+ Ar) = /(r) + A(r- V)/ + ]^\\r-Vff - \\\r-V)f + 0{\') (53) 
in Eq.(B2), the nullity of coefficients to the A and A^ powers implies, respectively, 

= :^ 111-212 [4/12 +(Rl2/l2)]dld2 (^4^) 

{BAh) 



1 
V 



6/i2 + ^(Ri2/i2) + ^(RL/i2: 



dld2 
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where we have adopted the operator notation given in (^Bj). Note that the sum of two 
zeros {Eq.(54a)}/8 + {Eq.(S46)}/4, 



V 



II - 21 



2/i2 + (Ri2/i2) + -(RL/: 



12 



dld2 



(i?5) 



with substitution /12 = F12 corresponds to the rhs of (|4y|) for the simplest N = 2 case. 
Let us introduce the functions gu and gu as follows 

TV 

Fugu = / d3 . . . diV n 



u,v — l 
(u<v) 

N 



„ N N 

Fugu = / d3...diV n 



{u<v) 



Both functions evidently possess the translational-invariance property. The sequence 
of operations {Eq.(55) for /12 = ^12^12} - {Eq.(54a) for /12 = Fi2(Ri2C/i2)}/4 + 
{Eq.(54a) for fu = F^guj/S leads to 
1 

V 







dld2 11 - 21 



2^12^12 + (Rl2i^l2)^ 



12 



g(Ri2-^i2)5'i2 + -Fugu 



(Rl2-Fi2)5'l2 + --^12(Rl25'l2) — --^12(Rl25'l2y 



We take advantage of the flexibility in choosing operator R (^61), and evaluate R^25'i2 first 
by applying R12 = (ri — r2) ■ Vi and then by applying R12 = (r2 — ri) ■ V2, with the 
result 



r 1 
Fi2{Rl29i2) = i"i2(Ri2^7i2) + i"i2(Ri2^i2) + Us . . . dN Y[ F^, -ET-ET 

u<v i,j=3 -^l*-^2j 
(Rli-^li)(R2j-^2i)} 



x{ [(r2 - n) • ViFi,] [(n - r,) ■ V2F2,- 
The substitution of (B9) into (B8) gives 







X 



i/ndzn^-ii-2i 

{2 + i-(R„F.,) + ^ 



[R12-F12) 



1 



N 



+ 2^ 

^ i=3 
1 1 



— {^liFiij + — (R2i-^2i) 



N r 



8F1 



[R12-F12) E 



+ 



12 

i,j=3 PliP2j 



i=3 



T^(Rli-^lj) + T^(R2i-^2j) 



(510) 



Rli-^li)(R2j-^2jj 



1 ^ 

N 



1 1 

- o E [(^2 - n) • ViFn] [(n - r,) ■ V2F2,] } 

Combining Eq.(BlO) for each pair of skeleton vertices with Eq.(^), the latter reduces to 
relation (48). 
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